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Abstract. We consider a quantum many-body system made of N interacting S—l/2 spins on a lattice, 
and develop a formalism which allows to extract, out of conventional magnetic observables, the quantum 
probabilities for any selected spin pair to be in maximally entangled or factorized two-spin states. This 
result is used in order to capture the meaning of entanglement properties in terms of magnetic behavior. In 
particular, we consider the concurrence between two spins and show how its expression extracts information 
on the presence of bipartite entanglement out of the probability distributions relative to specific sets of 
two-spin quantum states. We apply the above findings to the antiferromagnetic Heisenberg model in a 
uniform magnetic field, both on a chain and on a two-leg ladder. Using Quantum Monte Carlo simulations, 
we obtain the above probability distributions and the associated entanglement, discussing their evolution 
under application of the field. 

PACS. 03.67.Mn Entanglement production, characterization, and manipulation - 75.10.Jm Quantized 
spin models - 05.30.-d Quantum statistical mechanics 



1 Introduction 



Entanglement properties have recently entered the tool 
kit for studying magnetic systems, thanks to the insight 
they provide on aspects which are not directly accessi- 
ble through the analysis of standard magnetic observables 

PI21CI]0]|niinE] ' '^^ ie ana ly s ' s °t entanglement properties 
is particularly indicated whenever purely quantum effects 
come into play, as in the case of quantum phase tran- 
sitions. However, in order to gain a deeper insight into 
quantum criticality, as well as into other phenomena such 
as field-induced factorization and saturation, the con- 
nection between magnetic observables and entanglement 
estimators should be made clearer, a goal we aim at in 
this paper. On the other hand, most entanglement esti- 
mators, as defined for quantum magnetic systems, are ex- 
pressed in terms of magnetizations and spin correlation 
functions. It comes therefore natural to wonder where, in- 
side the standard magnetic observables, the information 
about entanglement is actually stored, and how entangle- 
ment estimators can extract it. Quite clearly, by posing 
this question, one does also address the problem of finding 
a possible experimental measure of entanglement, which 



is of crucial relevance in developing possible solid-state 
devices for quantum computation. 

In this context a privileged role is played by the con- 
currence C, which measures the entanglement of forma- 
tion between two q-bits by an expression which is valid not 
only for pure states but also for mixed ones 8,9. In the 
framework of interacting spin systems, exploiting differ- 
ent symmetries of such systems the concurrence has been 
related to spin-spin correlators and to magnetizations. |ll)l 
I11II12II13) However, C has not yet been given a general in- 
terpretation from the magnetic point of view, and a gen- 
uinely physical understanding of its expression is still elu- 
sive. 

Scope of this paper is therefore that of giving a simple 
physical interpretation of bipartite entanglement of forma- 
tion, building a direct connection between entanglement 
estimators and occupation probabilities of two-spin states 
in an interacting spin system. To this purpose we develop 
a general formalism for analyzing the spin configuration of 
the system, so as to directly relate it with the expression 
of the concurrence. The resulting equations are then used 
to read our data relative to the S — 1/2 antiferromagnetic 
Heisenberg model in a uniform magnetic field, both on a 
chain and on a two-leg ladder. The model is a cornerstone 
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in the study of magnetic systems, extensively investigated 
and quite understood in the zero-field case. When a uni- 
form magnetic field is applied the behavior of the system 
is enriched, gradually transforming its ground state and 
thermodynamic behavior. The analysis of entanglement 
properties in this model, and in particular that referring 
to the range of pairwise entanglement as field increases, 
sheds new light not only on the physical mechanism lead- 
ing to magnetic saturation in low-dimensional quantum 
systems, but also on the nature of some T = transitions 
observed in bosonic and fermionic systems, such as that 
of hard-core bosons with Coulomb interaction, and that 
described by the bond-charge extended Hubbard model, 
respectively. In the former case, the connection between 
magnetic and bosonic model is obtained by an exact map- 
ping that allows a straightforward generalization of our re- 
sults to the discussion of the phase diagram of the strongly 
interacting boson- Hubbard model |14| . In the more com- 
plex case of the the bond-charge extended Hubbard model, 
a direct connection between the Heisenberg antiferromag- 
nct in a field is not formally available, but a recent work 
by Anfossi et al.]U^ has shown that some of the T = 
transitions observed in the system are characterized by 
long-ranged pairwise entanglement of the same type we 
observe in our magnetic model at saturation. 

Our data result from stochastic series expansion (SSE) 
quantum Monte Carlo simulations based on the directed- 
loop algorithm|lf>|. The calculations were carried on a 
chain with size L = 64 and on a L x 2 ladder with L = 40. 
In order to capture the ground-state behavior we have 
considered inverse temperatures (3 = 2L. 

In SccFJlwc define the magnetic observables we refer to, 
and develop the formalism which allows us to write them 
in terms of probabilities for two spins to be in specific 
states, both at zero and at finite temperature. In Sec. 
we show how concurrence extracts, out of the above proba- 
bilities, the specific information on bipartite entanglement 
of formation. In Sees, rj] and we present our SSE data for 
the antiferromagnetic Heisenberg model on a chain and 
on a square ladder respectively, and read them in light of 
the discussion of Sees. El and Conclusions are drawn in 
Sec. El 



2 From magnetic observables to spin 
configurations 

We study a magnetic system made of N spins 5 = 1/2 
sitting on a lattice. Each spin is described by a quantum 
operator Si, with [Sf,S!^\ = iSi m e a f3jS^ , I and m being 
the site-indexes. 

The magnetic observables we consider are the local 
magnetization along the quantization axis: 



(i) 



and the correlation functions between two spins sitting on 
sites I and m: 



The averages ( • ) represent expectation values over the 
ground state for T = 0, and thermodynamic averages for 
T > 0. 

We now show that the above single-spin and two-spin 
quantities provide a direct information on the specific quan- 
tum state of any two spins of the system. Let us consider 
the T = case first: For a lighter notation we drop site- 
indexes, allowing their appearance whenever needed. After 
selecting two spins, sitting on sites I and m, any pure state 
of the system may be written as 



(3) 



ren 



where S is an orthonormal basis for the 4-dimensional 
Hilbert space of the selected spin pair, while 1Z is an or- 
thonormal basis for the 2 Ar " 2 -dimensional Hilbert space of 
the rest of the system. Moreover, in order to simplify the 
notation, we understand products of kets relative to (oper- 
ators acting on) different spins as tensor products, mean- 
while dropping the corresponding symbol <8. The quantum 
probability for the spin pair to be in the state \v), being 
the system in the pure state \&), is p v = J2 r \c„ r \ 2 , and 
the normalization condition l^^) = 1 implies ~^2 v Pv = 1. 
We consider three particular bases for the spin pair: 



with 



S 1 ={\u l ),\u I ),\u ]l[ ),\u r ,)} 
<S 2 ={|ei), |e 2 ), |e 3 ), |e 4 » , 
<S3={|u,), K), |e 3 ), |e 4 }} , 



K) = I T>i| T)™ , K) = I I) 



(4) 
(5) 
(6) 



mi: 



m i | "Tv/ 



D/IT) 



m ) 
m i 



lei) = 7! (K) + K» , |e 2 ) = 75 (K) - K» , 
|e 3 ) = 7= (K) + \u w )) , |e 4 ) = 7= (K) - \u w )) , (7) 

where | T)i,ro(| i)l,m) are eigenstates of 5f m with eigen- 
value +§(— For the coefficients entering Eq. ©, and 
for each state r, the following relations hold 

k ( C ir + c nr) > C 2r = "75 ( C ir ~ c .r) » ( 8 ) 



C 3r ( C 'r C wr) ' C 4 

meaning also 



V2 

71' 



O , (9) 



kuf + k r | 2 = M a + Krl 2 > 



(10) 

\c ir \ 2 - \c 2r \ 2 = |c Ir c nr |cos(^ r - ^) , (11) 
|c 3r | 2 + |c 4r | 2 = |c mr | 2 + | Clvr | 2 , (12) 

|c 3 r| 2 - |c 4r | 2 = |c mr C IVr |cOs(^ - If?) , (13) 

where c vr = \c^ r \e l{p ^ . 

According to the usual nomenclature Si and S2 are 
the standard and Bell bases, respectively, while S3 is here 
called the mixed basis. Such bases are characterized by the 
fact that states corresponding to parallel and antiparallel 
spins do not mix with each other. It therefore makes sense 
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to refer to |M r ), |w e ), |ei), and \e2) as parallel states, and to 
l M m)j I u iv)j l e 3)j ancl l e 4) as antiparallel states. The proba- 
bilities specifically related with the elements of 5i will be 
hereafter indicated by P^PhtPw and P N while pi,p2,P3, 
and pi will be used for those relative to the elements of 
1S2. From the normalization conditions 

Pi +Pe +Pm +Pw = 1 ( 14 ) 
Pi +P2 +P3 +P4 = 1 (15) 
Pi + P k +P3 +P4 = 1 , (16) 

or equivalently from Eqs. (jlOJI and (|12Jl . follows j3 : +p n = 
Pi + _P2, and p ra + p w = p 3 + P4, representing the prob- 
ability for the two spins to be parallel and antiparallel, 
respectively. We do also notice that the elements of Si are 
factorized states, while those of 1S2 are maximally entan- 
gled ones. 

The above description is easily translated in terms of 
the two-site reduced density matrix 

P = £<r|W|r) = ]>>)(A|5> rC ; r , (n) 

r i>\ r 

whose diagonal elements are the probabilities for the ele- 
ments of the basis chosen for writing p. The normalization 
conditions Eqs. (|14I16() translate into Tr (p) = 1. 

Thanks to the above parametrization, the magnetic 
observables and J2Jl are directly connected to the prob- 
abilities of the two spins being in one of the states J7J. In 
fact it is 
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= (*\s+s- + SfS+\ (| e3 )5> 3r |r} + | e4 )5> 4r |r> ) 



<#i ^ 3 )E c 3.i r )-i e 4)E c -i r )) 



(pa -Pa), (18) 
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and similarly 



to the notation of the previous section, it is easily shown 
that 



2{g xa 

9 ZZ = 
M z = 



g yv ) = (pi-P2) 



(19) 



3(p.+a)-s = 5 (pi +P2)-h ( 2 o) = K c ' ~ - ( c ' ~ c :)l = 2 l ( 



CttT 



(30) 



|(M,« + M*) = (p I - A ) , 



(21) 



where all 5i are suitable to calculate g zz , while (g xx ±g vy ) 
and M z specifically require 52 and 1S3, respectively. After 
Eqs. one finds 



Pi = 


1 

4 


+ g zz + M z 


! 


(22) 


Pn = 


1 

4 


+ 9 ZZ - M z 


5 


(23) 


Pi = 


1 

4 


+ g xx - g vv 


+ g zz , 


(24) 


P2 = 


1 

4 


- g xx + gvv 


+ g zz , 


(25) 


P3 = 


1 

4 


+ g xx + 9 VV 


-g zz , 


(26) 


Pi = 


1 

4 


- g xx - g vv 


-g zz ■ 


(27) 



It is to be noticed that the probabilities relative to the 
Bell states do not depend on the magnetization. 

In the the finite temperature case, the generalization is 
straightforwardly obtained by writing each of the Hamil- 
tonian eigenstates, numbered by the index n, in the form 
©, so that 



p(r) = 5>M£ e "* ,/T £' 



(28) 



In terms of probabilities the above expression simply means 
that the purely quantum p^ shall be replaced by the quan- 
tum statistical probabilities 



p^-E^Ei 



(29) 



Therefore, apart from the further complication of the for- 
malism, the discussion developed for pure states stays sub- 
stantially unchanged when T > 0. 

Equations Q22 [l -Q27 jl show that magnetic observables 
allow a certain insight into the spin configuration of the 
system, as they give, when properly combined, the prob- 
abilities for any selected spin pair to be in some spe- 
cific quantum state. However, the mere knowledge of such 
probabilities is not sufficient to appreciate the quantum 
character of the global state, and more specifically to quan- 
tify its entanglement properties. 



3 From spin configurations to entanglement 
properties 

We here analyze the entanglement of formation 17,8,9 
between two spins, quantified by the concurrence C. In the 
simplest case of two isolated spins in the pure state \<f>) the 
concurrence may be written as C = | J2i a i I > where on are 
the coefficients entering the decomposition of \4>) upon the 
magic basis {|ei), i|e2), i\e-s), However, if one refers 



where Eqs. JSJ-lO have been used, with index r obviously 
suppressed. The above expression shows that C extracts 
the information about the entanglement between the two 
spins by combining probabilities and phases relative to 
specific two-spin state. 

In fact, one should notice that a finite probability for 
two spins to be in a maximally entangled state does not 
guarantee per se the existence of entanglement between 
them, since this probability may be finite even if the two 
spins are in a separable state. |T%] In a system with de- 
caying correlations, at infinite separation all probabilities 
associated to Bell states attain the value of 1/4, but this 
of course tells nothing about the entanglement between 
them, which is clearly vanishing. It is therefore expected 
that differences between such probabilities, rather than 
the probabilities themselves give insight in the presence 
or absence of entanglement. 

When the many-body case is tackled, the mixed-state 
concurrence of the selected spin pair has an involved defi- 
nition in terms of the reduced two-spin density matrix. 9 
However, possible symmetries of the Hamiltonian TL greatly 
simplify the problem to the extent that C results a simple 
function of the probabilities l|22|l - (|27JI oruv - We here as- 
sume that Ti. is real, has parity symmetry (meaning that 
either Ti. leaves the z component of the total magnetic 
moment unchanged, or changes it in steps of 2), and is 
further characterized by translational and site-inversion 
invariance. The two latter properties implies as de- 
fined in Eq. (TjQl to coincide with the uniform magnetiza- 
tion M z = ^2i(S z )/N, and the probabilities p m = p w , 
respectively. 

Under these assumptions, the concurrence for a given 
spin pair is|13j 



C\ r) =2max{0,C ( V ) ,C ( '' } 



(r) 



a 



\„xx , yy I 

\9( r ) + g {r) \ 

yy 



9'i 



M 2 



(r) - l.9(r) g { r) 



g{; 



(31) 

(32) 
(33) 



where r is the distance in lattice units between the two 
selected spins. Despite being simple combinations of mag- 
netic observables, the physical content of the above ex- 
pressions is not straightforward. However, by using the 
expression found in Section |2 one can write Eqs. i|32|) and 
(I33|l in terms of the probabilities for the two spins to be in 
maximally entangled or factorized states, thus finding, in 
some sense, an expression which is analogous to Eq. (|30|l 
for the case of mixed states. In fact, from Eqs. (|18fl - l|19|) . 
it follows 



2(7'= |ps-p 4 | . 
2C" = \pi - p 2 \ - (1 - pi - P2) 

= |Pl -P2| - VPmPiv > 



(34) 
(35) 
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where we have used p m =p w and hence P3 + Pi = 2p m = 
2^/p m p w . The expression for C may be written in the 
particularly simple form 

2C" = 2max{pi,p 2 }-l , (36) 

telling us that, in order for C" to be positive, it must be 
either p\ > 1/2 or p2 > 1/2. This means that one of the 
two parallel Bell states needs to saturate at least half of 
the probability, which implies that it is by far the state 
where the spin pair is most likely to be found. 

Despite the apparently similar structure of Eqs. 134f) 
and understanding C is more involved, due to the 
fact that ^fPiP^ cannot be further simplified unless p : = 
. The marked difference between C and C" reflects the 
different mechanism through which parallel and antiparal- 
lel entanglement is generated when time reversal symme- 
try is broken, meaning p 1 ^ p n and hence M z 7^ 0. In fact, 
in the zero magnetization case, it is p E = p 1 = +p2)/2 
and hence 

2C" = 2max{p 3 ,p 4 }-l , (37) 

which is fully analogous to Eq. (|36|l . so that the above 
analysis can be repeated by simply replacing p\ and P2 
with P3 and P4 . 

For M z ^ 0, the structure of Eq. (|37|) is somehow kept 
by introducing the quantity 

A^^-^f, (38) 

so that 

2C' = 2max{p3,p 4 } - (1 - A 2 ) , (39) 

meaning that the presence of a magnetic field favors bi- 
partite entanglement associated to antiparallel Bell states, 
\e%) and \ei). In fact, when time reversal symmetry is bro- 
ken the concurrence can be finite even if p%, p^ < 1/2. 

From Eqs. I|36|l and (|39|l one can conclude that, de- 
pending on C being finite due to C or C", the entangle- 
ment of formation originates from finite probabilities for 
the two selected spins to be parallel or antiparallel, re- 
spectively. In this sense we will speak about parallel and 
antiparallel entanglement. 

Moreover, from Eqs. (|34|I - H35|I we notice that, in order 
for parallel (antiparallel) entanglement to be present in 
the system, the probabilities for the two parallel (antipar- 
allel) Bell states must be not only finite but also different 
from each other. Thus, the Bell states |ei) and |e 2 ) (|e3) 
and I £4)) result mutually exclusive in the formation of 
entanglement between two spins in the system, the latter 
being present only if one of the Bell state is more probable 
than the others. The case pi — p 2 — 1/2 (p 3 = p 4 = 1/2) 
corresponds in turn to an incoherent mixture of \e\) and 
|e 2 ) (|e 3 ) and |e 4 )). 

In fact, the occurrence of the differences \pi — p 2 | and 
\P3 — P4 1 is intriguing. Let us comment on \pi — as the 
same kind of analysis holds for \p% —p 4 \ . In the general case 
the difference p\ —P2 can vanish because of genuine many- 
body effects which are not directly readable in terms of 2- 
spin entangled or separable states. It is easier to interpret 
Eq. (|35(l [Eq. I)34[l]. if one restricts the possibilities to the 
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Fig. 1. Magnetization and correlators versus the magnetic field 
h for the chain Eq. 1401 . The dashed lines mark the value of 
the field where p n = (see text). 

case in which the two spins are not entangled with the 
rest of the system. By using Eq. (|1 1|) . one can select two 
particular situations all leading to p\ = P2 : 

(i) c lr or c Er vanishes VT, meaning that \&) does not 
contain states where the two selected spins are parallel 
and entangled; 

(ii) for each r such that both \c lr \ and |c nr | are non- 
zero, it is ip^ — (p^ = 7r/2. Thus, whatever the antiparallel 
components are, the parallel terms of \i r ) appear in the 
form (a|ei) + a*|e 2 )). 

The above analysis suggests the first term in C" (C) 
to distill, out of all possible parallel (antiparallel) spin 
configurations, those which are specifically related with 
entangled parallel (antiparallel) states. These characteris- 
tics reinforce the meaning of what we have called parallel 
and antiparallel entanglement. 

4 Chain 

We consider the isotropic Heisenberg antiferromagnetic 
chain in a uniform magnetic field, described by 

j = Y i S i -S i+1 -hS? , (40) 

i 

where the exchange integral J is positive, and the reduced 
magnetic field h=g^H / J is assumed uniform. 

This model is characterized by the rotational symme- 
try on the xy plane, as well as by the existence of a satu- 
ration field h s = 2, such that for h > h s the ground state 
is the factorized ferromagnetic one, with all spins aligned 
along the field direction. Moreover, Eq. (|40|l has all the 
necessary symmetries for Eqs. <|31[1 - <|33|) to hold. 

Due to the rotational symmetry on the xy plane, it is 
gxx _ gyy^ meanm g p 1 = p 2 = I _|_ < 1/2, according to 

Eqs. (|24() and l|25() . and hence null parallel entanglement 
(C" < 0) between any two spins along the chain, no mat- 
ter the field, the temperature, and the distance between 
them. 

In Fig. n we show the T = correlation functions for 
nearest neighboring (n.n.) and next-nearest neighboring 
(n.n.n.) spins, together with the uniform magnetization, 



Andrea Fubini et al.: Reading entanglement in terms of spin configurations in quantum magnets 



7 



0.75 



0.25 




0.75 




0.25 



Fig. 2. Concurrence and related probabilities of the mixed Fig. 3. Concurrence and related probabilities of the mixed 
basis states S3 [Eq. @] for n.n. sites of the chain Eq. (14011 basis states S3 [Eq. JSJ] for n.n.n. sites of the chain Eq. (1401 . 



as the field is varied. Beyond the overall regular behavior, 
we notice that there exists a value of the magnetic field 
where one simultaneously observes gf^. = and M z = 1/4 
(indicated by the dashed lines). According to Eqs. 1|2U[I 
and (|21|l this implies null probability p n for adjacent spins 
to be parallel in the direction opposite to the field. This 
means that the ground-state configuration is a superposi- 
tion of spin configurations entirely made of stable clusters 
of spins parallel to the field separated by Neel-like strings. 

In Fig. [21 we show the probabilities for n.n. spins to 
be in the states of the mixed basis, together with the n.n. 
concurrence: The value of the n.n. concurrence for h = is 
in agreement with the exact resut in the thermodynamic 
limit .[H] In presence of an external magnetic field, Cm 
is found positive V/i, meaning that, no matter the value 
of the field, the probabilities P3 and p^ for adjacent spins 
are always different from each other. The probabilities for 
the triplet states \e$), |u n ), and are equal for h = 
and depart from each other when the field is switched on. 
The singlet |e4) evidently dominates the ground state up 
to a field which roughly corresponds to the value where 
p s vanishes. 

As for the concurrence, despite the ground-state struc- 
ture evidently changes as the field increases, Cm stays 
substantially constant up to a large value of the field, 
mainly due to the fact that not only p± but also ps de- 
creases with the field. This behavior mimics the one oc- 
curring in a spin dimer, whose ground state is the singlet 
state \ei) up to h = 1 where, after a level crossing, 
becomes energetically favored. However, in a spin chain, 
many-body effects smear the sharp behavior of the dimer 
due to the level crossing. We do also notice that Cm starts 
to decrease as soon as the total probability for parallel 
spins (p 1 + p n ) gets larger than that for antiparallcl spins 
{P3 + Pi)- The further reduction of C(i) is mainly driven 
by p I starting to rapidly increase. 

In the same field region where a substantial change 
in the n.n. configuration occurs, the n.n.n. concurrence 
C( 2 ) switches on. This is seen in Fig. [31 where the proba- 
bilities for n.n.n. spins arc shown together with the corre- 
sponding concurrence. In fact, when considering the n.n.n. 
quantities, we notice that both gf£ and gf^, have a non- 
monotonic behavior, displaying a maximum and a mini- 
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Fig. 4. From the upper to the lower curves: concurrences from 
the nearest- up to the 5 th -neighbors of the chain Eq. 14U1 . The 
inset shows the divergence of the range of the concurrence as 
h — > h s , the line shows the (h — /i s )~ 1,/2 behavior. 



mum, respectively, in the field region where C( 2 ) gets pos- 
itive (as from the comparison between Fig. [T] and 0. 

Regarding the probabilities, one finds that, although 
the most likely state is always p^ is surprisingly large, 
and almost equal to p T , as far as h < 1. Moreover, both p$ 
and pi have a non monotonic behavior and increase with 
h up to the field where we simultaneously observe g^ and 
<?H attaining their extreme values, p 1 exceeding 1/2, p^ 
getting larger than p n , and C( 2 ) switching on. 

As observed in the n.n. case, when p n for n.n.n. spins 
vanishes C/3) switches on. Let us further comment upon 
C(i), C(2), and Cr^y Given the fact that only antiparallcl 
entanglement may exist in this chain, it is not surprising 
that Cm > and Cn) — at low fields, as n.n. spins 
belong to different sublattices, while n.n.n. spins belong 
to the same sublattice. However, the fact that C(2) be- 
comes finite indicates a ground-state evolution from the 
Neel-like to the ferromagnetic state such that the system 
enters a region where quantum fluctuations increase the 
total probability for spins belonging to the same sublattice 
to be antiparallel and entangled. The opposite effect is un- 
derstood when C(3) is considered: in order to keep C(3) = 
almost up to the saturation field, quantum fluctuations 
must reduce the total probability for spins belonging to 
different sublattices to be antiparallel and entangled. 

The above comments upon C( 2 ) and C/3) may be gener- 
alized to C(„) with even and odd n, respectively. In Fig.0] 
we in fact show C n up to n — 5. The concurrence for in- 
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creasing distance between the two spins gets finite for a 
big enough field resembling the phenomenology of finite 0.6 
spin clusters. £Q Moreover, combining the exact results of 
Refs.^Hl and [23j, we find that the range of the concur- 
rence for the model l|4(J|) , namely the distance R such that 0.4 
C/ r ) vanishes for r > R, is 



R = 



V^(2 + 4M z )(i-M z )V2 



(41) 



with the constant p = 0.924.... When h — > h a , it is M z ~ 

| - ^rVh s - h and 6 ~ 2 - ^V^s - h, and the range 
of the concurrence is seen to diverge according to R ~ 

^ff-{h s — h)^ 1 / 2 . In other terms, approaching the satura- 
tion field, all C(„) become finite of order 0(1/ N), con- 
sistently with the occurrence of a |Wjv) state|2]. F° r 
such state the entanglement is maximally bipartite in the 
sense of the Coffman-Kundu-Wootters inequality. |2 111221 
I23| This scenario is consistent with our numerical data. 
As shown in FigQ] up to n = 5, for any C(„) it exists a 
field h n > h n -\ such that C(„) is positive for h £ [h n ,2), 
with h n — ► 2 for n — ► oo. The divergence of the range of 
the concurrence for /i — + h s is shown in the inset of Fig. 
Although the correct power-law behavior shows up, the 
precision of the numerical data is not sufficient to get the 
correct multiplicative constant. In fact, the above expres- 
sion l|41(l is derived from asymptotic exact results, valid 
only for r> 1, when Cr r ) becomes too small to resolve it 
numerically. 

The formalism introduced in the previous sections works 
also in the finite temperature case, where it describes the 
effects of thermal fluctuations on quantum coherence. In 
Fig. the temperature dependence of probabilities and 
concurrences, for h = 1.8, shows how thermal fluctua- 
tions progressively drive the system towards an incoher- 
ent mixture of states. Increasing T the concurrences (right 
panel) are progressively suppressed and above k B T ~ 0.8 J 
also the n.n. concurrence vanishes. At higher temperatures 
none of the spin pairs in the system is entangled and quan- 
tum coherence is lost. The temperature behavior of the 
probabilities (left panel) is non monotonic, signaling the 
relative weight of the different states in the energy spec- 
trum of the system. Eventually, at high T all the proba- 
bilities tends to the asymptotic value p v — 1/4. 



5 Two- Leg Ladder 

The above picture further enriches when considering the 
two-leg isotropic ladder, described by 

7 = E E ' S i+h*-hSl a ) +7^,0 • S M , (42) 

t a=0,l 

where the index i runs on both the right (a = 0) and 
left (a = 1) leg. The first term is the Heisenberg Hamil- 
tonian l|4U|l for the right and left legs, while the last term 
describes the exchange interaction between spins of the 
same rung, whose relative weight is 7. 
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Fig. 5. Left panel: n.n. probabilities versus temperature at 
h — 1.8. Right panel: n.n. and n.n.n. concurrences versus tem- 
perature at h — 1.8. 



The model (|42|) is known to describe cuprate compounds 
like SrCu2C>3 and it has been extensively studied for zero|2"H 
12*5] and finite field [23 • The system shows a gap A in the 
excitation spectrum that can be interpreted essentially as 
due to the energy cost for producing a triplet excitation on 
a rung [23 ■ The system reaches full polarization[2E] , with 
all spins aligned along the field direction, for h > 7 + 2. 

In the following we will specifically consider the isotropic 
case 7 = 1, which is characterized by a gap A ~ 0.5.7. |24| 
and by a saturation field h s — 3. As in the chain case, due 
to the rotational invariance on the xy plane, parallel en- 
tanglement cannot develop in the isotropic ladder. On the 
other hand, antiparallel bipartite entanglement can here 
develop between spins belonging to the same leg, or to the 
same rung, or to a different rung and leg. Two-spin quan- 
tities will be hereafter pinpointed by the two-component 
vector (n, r a ) joining the two selected spins, the first com- 
ponent referring to the direction of the legs, and the second 
one to that of the rungs. The indexes (01), (10), (11), (20) 
will therefore indicate n.n. spins on the same rung, n.n. 
along one leg, n.n.n. on adjacent rungs, and n.n.n. along 
the same leg, respectively. 

Our SSE data in Fig. for the uniform magnetiza- 
tion and the n.n. correlation functions 9?^y 9nm connrm 
the description given in the previous paragraph: Before 
the Zeeman interaction fills the energy gap at the critical 
value h c ~ 0.5, the ground-state configuration is frozen 
and characterized by the singlet \e4) being by far the most 
likely state for each rung. 

The use of the formalism developed in Sec. gives a 
direct information on the physics of the system: In Fig. 
we see that the singlet probabilities P4 relative to n.n. spins 
on a rung and along one leg, as functions of the field, share 
a similar behavior everywhere but at the critical field h c , 
where for n.n. spins sitting on the same rung shows 
up a kink that is not present in the singlet probability 
along the leg. This qualitatively different behavior clearly 
reflects the nature of the energy gap that closes at h c . 
The sharp decrease of in favor of p l on the rung just 
above h c testifies that, even in the case, here considered, 
of equal exchange interaction along the legs and on the 
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Fig. 6. Magnetization and correlators versus the magnetic field 
h for the ladder Eq. 
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Fig. 7. Probabilities relative to the mixed basis versus the 
magnetic field h for the ladder Eq. 14211 : p : (0), p^ (□), P3 
(A), and P4 (V). Open (full) symbols are for n.n. spins along 
the same leg (on the same rung). The inset zooms in on the 
behavior of p r , p E , and p3 near the critical field h c . 



rungs (7 — 1), the first excitations in the energy spectrum 
of the ladder are triplet excitations on the rungs. 

When the Zeeman energy becomes larger than the gap, 
for h > h c , the ground state starts to evolve with the field, 
whose immediate effect is that of pushing the quantities 
relative to spins on the rungs and along the legs towards 
each other: In fact, for h > 1 n.n. spins along the legs and 
on the rungs substantially share the same behavior. As for 
the probabilities, we see that p n and P3 keep being equiv- 
alent, no matter the value of the field, and slowly vanish 
as saturation is reached. On the contrary the probability 
for n.n. spins to be in increases at the expense of the 
probability relative to the singlet state until, for h ~ 1.8, 
the two probabilities cross each other. 



' : same leg 

0.3 '- 




Fig. 8. From the upper to the lower curves: concurrences rela- 
tive to spins belonging to the same (upper panel) and different 
legs (lower panel) versus field for the isotropic ladder up to 
r = 4. 

Finally, we apply the formalism of Sec. 01 to extract 
features of the ground state from the concurrences. Fig. El 
shows the concurrences Cr rt ra \ up to the distance r = 
T"i+r a = 4 for spins sitting on the same (upper panel) and 
on different legs. The bipartite antiparallcl entanglement 
between two spins sitting at a given distance r is in general 
larger on different legs, even beyond the n.n. case. 

As expected, the field, after closing the gap, pushes 
C(oi) and Cnu) towards each other. Quite unexpectedly, 
however, this evolution includes a region where n.n. con- 
currence along the leg, C(i ), increases. It is interesting to 
notice that C(\\\ switches on at h ~ 1.8, where p^ and 
p I for n.n. spins are seen to cross each other in Fig. |6] 
signaling the crossover from an antiferromagnetic to a 
ferromagnetic-like configuration of the n.n. spins. 



6 Conclusions 

In this paper we developed a simple and effective formal- 
ism that allows to reconstruct the probability for two spins 
of a multi-spin system to be in a given quantum state, once 
the collective state of the system is given. Remarkably, 
such probabilities are found to be simple combination of 
standard magnetic observables, Eqs. (|22I27(1 . Within such 
formalism it is very natural to understand how concur- 
rence quantifies the amount of entanglement between two 
spins by comparing the probabilities for those spins to be 
in different Bell states. In particular the expression for the 
concurrence clearly separates the case of parallel [Eq. 1)36(1 ] 
and antiparallel [Eq. ffify ] spins, leading to the introduc- 
tion of the concept of parallel and antiparallel entangle- 
ment. 

The knowledge of the probability distribution for a 
given set of two-spin states can be a useful tool to study 
quantum phases dominated by the formation of particu- 
lar local two-spin states and to investigate the transitions 
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given by the alternation of such states. Within this class 
of phenomena we can cite the occurrence of short-range 
valence-bond states in low-dimensional quantum antifer- 
romagnets and the transition from a dimer-singlet 
phase to long-range order in systems of weakly coupled 
dimcrs under application of a field or by tuning of the 
inter-dimer coupling |28| . 
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